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We compute the dependence of the tunneling current in a double point contact in the k = 3 
Read-Rezayi state (which is conjectured to describe an incompressible quantum hall fluid at filling 
fraction v = 12/5) on voltage, separation between the two contacts, and temperature. Using the 
tunneling hamiltonian of cond-mat/0607431 , we show that the effect of quasiholes in the bulk region 
between the two contacts is simply an overall constant multiplying the interference term. This is the 
same effect as found for the differential conductivity in cond-mat/0601242 the difference is that we 
do an actual edge theory calculation and compute the full current-voltage curve at weak tunneling. 



I. INTRODUCTION 



Certain quantum systems in 2 dimensions can support quasiparticles with anyonic and non-abelian statistics. 
These are experimentally realized in fractional quantum hall (FQH) systems. Specifically, the v = 5/2 FQH state 
is conjectured to be in the universality class of the Moore-Read (MR) Pfaffian wavefunction, and it is believed that 
there is a FQH state at v — 12/5 described by the k = 3, M =1 Read-Rezayi (RR) state. Both are non-abelian, and 
the latter supports universal quantum computation, making its study especially important. 

In this paper we consider a double point contact geometry in an assumed v — 12/5 RR state, and study tunneling 
in the edge theory. A similar setup is considered in [l[. The edge theory consists of a charged boson, and a neutral 
sector described by the Z3 parafermion CFT. The latter contains spin fields <7i which are not local with respect to the 
other fields in the CFT, and it is precisely these that enter the hamiltonian describing tunneling of quasiholes, which 
is the most relevant tunneling process at low temperatures. Unfortunately it is difficult to interpret such a tunneling 
Hamiltonian, because in a chiral theory, correlation functions of spin fields are not well defined. Instead, they could 
be any linear combination of a finite number of conformal blocks. In fact, it turns out that the choice of such linear 
combination is intimately related to the topology of the corresponding tunneling process occurring in the bulk; Q 
give a prescription for determining the correct choice of linear combination, in the case of a single point contact. We 
adapt their method to the case of a double point contact and use it to compute the tunneling current to first order 
in perturbation theory. In contrast to the abelian case^, we find that the interference term in the tunneling current 
is suppressed depending on the topological charge contained in the bulk region between the two point contacts. A 
crucial point is that a bulk quasihole affecting the tunneling current must be prepared and annihilated on the edge 
in the far past and future respectively, and the prescription of Q must be applied to this entire process. After one 
divides out the universal factor associated with preparing and annihilating the bulk quasihole, one obtains correlation 
functions of the original spin fields which are finite and correctly incorporate the effect of the bulk quasihole. 

II. SETUP 

We assume that there exists an incompressible quantum hall fluid at v = 12/5, and that it consists of a filled 
lowest Landau level (of both spins) and a first Landau level described by the fc = 3, M = 1 Read-Rezayi state, which 
has filling factor 2/5. The Read-Rezayi wave function is a correlator of operators in the conformal field theory of a 
charged boson and Z3 parafermions. The parafermions contain 6 fields: 1, a%, tpi, e, and their fusion rules are 
determined by o\(72 = 1 + e and ipiip2 = 1 0- We will need the electron and quasihole operators; they are as follows 3 : 

$ cl . = ^eVI^ (1) 
$ q .h. = <r ie vM<= (2) 

We have q c \. = 1,A^ )1 = |jA e i. = §, and g q .h. = g,A CT1 = yg,A q .h. = ^ . The chiral edge theory also consists of 
Z3 parafermions and a charged mode, but they propagate with different velocities, v n and v c . Generally, one expects 
v c > v n . In writing down the edge theory description of the tunneling hamiltonian, we a priori need to consider all 
sorts of tunneling processes, but a quick renormalization group argument 2 - shows that the most relevant term at low 
temperatures corresponds to an e/5-charge quasihole - the tunneling hamiltonian simply destroys such a quasihole 
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on one side of the point contact and creates it on the other. Here we assume we are at temperatures higher than 
those at which the crossover physics of 0, [1], etc. becomes important, so that the tunneling is weak and can be 
treated perturbatively. This is an important point, because when we make reference to zero temperature results from 
now on we actually mean temperature small compared to eV, where V is the bias voltage, and to ft^osc an d ^ w "sc- 
Here w„ sc = 2a /v c and w" sc = 2a /v n (a being the distance between the two contacts) but still large compared to the 
crossover temperature. Note that we have 4 different energy scales, set by the temperature, bias voltage, wj; sc and 
w™ sc . The last two cannot be varied independently, but only through the variation of the separation a. The two cases 
we consider in this paper are zero temperature, where effectively we have one dimensionless parameter we can vary 
(say eV / (hu>o SC )) , and temperature high compared to huj osc , where things simplify as well. 
The tunneling hamiltonian is 



H t (t) = yY j e^ t V j {t) + h. c. (3) 



where Tj is the tunneling amplitude at the jth point contact (j — 1,2), and ujj = |Jr is the Josephson frequency 
corresponding to a charge e/5 quasihole. The tunneling operator Vj is defined as 



Vj(t) = ^(^, I ,^ 2 fe 1 ,,t)e*^ ( ^" t) e-*^ (l -' t) (4) 

where Xjj and Xj tU are the spatial coordinates of the lower and upper sides of the jth point contact. Note that because 

the hall bar is large but finite, the upper and lower points arc actually on the same edge, just very far away from each 
other. From now on it will be convenient to use the parametrization 

x j,i — x j (5) 

Xj, u = D — Xj (6) 

where D is the distance around the hall bar, which we will eventually take to infinity. 



III. TUNNELING CURRENT IN PERTURBATION THEORY 

Our discussion here closely parallels , except that we have the added complication of dealing with the parafermions 
and spin field operators. As mentioned before, there is ambiguity in treating the <Ji as operators, and their correlators 
are not well defined. Indeed, to resolve the ambiguity we will need to add extra Oi operators, which create and 
annihilate the bulk quasiholes in the far past and future. As we shall see, taking these to infinity essentially decouples 
them from the rest of the correlator, and their only influence is topological. In fact, together with making the hall 
bar really large (that is, taking D to infinity), the correlator of interest factors as a product of 2-point functions, 
which is determined by conformal invariance up to a multiplicative constant. All of the effects mentioned above are 
incorporated in that constant, denoted by Cjk below. In this section we calculate the tunneling current leaving Cjk 
undetermined. In the next section we determine Cjk- 

We begin by writing down the appropriate operators. The operator for the current across the point contact is 
the commutator of the total charge and the tunneling hamiltonian; because only the charge mode is involved the 
calculation is easy and we get 



j® = J n^Tje-^Vj^-h. c.j (7) 
To lowest order in perturbation theory, the expectation value of the tunneling current is 



We have 



(?(*)>=-< f dt'{0\\j(t),H t (t')]\Q). 



(8) 
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W)H t {t')) = ^^((iy-^^Ci) -iy e ^V/(i)) (r fee -^ + r fc v^ t V fc t (i'))) 

= | E( r J r ^ e ^ J(t ^' ) ^W^ t ( < ')) -r/r^'"''^/^^))) (9) 



because the only terms that contribute to the vacuum expectation value are ones that conserve the total charge on 
each side of the hall bar. Now, 

a x (D - x kt t')a 2 (x k , ^-^(^O^^d^,^ (10) 

= C jk (cri(xj,t)a2(xk,t'))(a2(D - Xj,i)(j\{D - x k ,t')) 

(e^^ ( ^'* ) e "^ 0c(afc '*' ) )(e"^ 0c(£ '" :,;j ' t) e^s </ ' <:(I) ~ ;Efe ' t ' ) ) (11) 

= C jk [6 + i(v n (t - + {xj - x k ))r 2/15 {6 + i(v n {t - t') - (x 3 - x k ))}- 2/15 

[S + i{v c (t - t') + ( Xj - x k ))r 1/15 [S + i(v c (t - t') - ( Xj - x k ))]- 1/15 . (12) 

Here Cj k is the constant mentioned above, which is determined by choice of conformal block (its phase is more subtle 
and convention dependent), and <5 is an infinitesimal, specifying the choice of branch cut. We will obtain all the 
information about Cj k we need for our purposes in the next section. For now, define 

pa, x) = [s + i(v n t + x )r 2/i5 [s + i{v n t - x )r 2/15 [s + i{v c t + x)r i/15 [s + i{ Vc t - x)r i/i5 . (13) 

so that 

(V^V^it')) =C jk P(xj-x k ,t-t'). (14) 

In a similar manner we obtain 

(V^(t)V k (t')) = C' ]k P{x -x k ,t- t') = C' jk P(x k ~ Xj ,-(t' - t)) (15) 

Because Vj'(t) and V k (t) commute for spacelike separation, we have Cj k = C' k j (for j — k the argument is a little 
more subtle but the conclusion holds - basically you have to infinitesimally spacelike separate the two operators and 
then argue using continuity). We therefore have 

(\j(t), H t (t')}) = ^V"(T j r k *e-^ t -^C jk [P(x j -x k ,t-t')-P(x j -x k ,-(t-t'))] 
5 j,k 

-rfrte^rt-^CkjiPixk -Xj,t- t') - P{x k - x,, -(t - t'))}) (16) 
Doing the integral in © we obtain 

OX*)) = ^E^W + CkjTfT^PiujtXj - x k ) - P{-uj, Xj - x k )\ (17) 

j-k 

where P(loj,x) is the Fourier transform of P(t,x). For x = we can do the integral explicitly: 

P( W> 0) = ^- 4/1 V- 2/15 f ^ykr 3/5 6(c) (18) 
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When x 7^ 0, we can write down an alternate integral representation for P(t,x) as follows: in the integral for the 
Fourier transform, the denominator is a product of 4 terms, with t (the variable being integrated over) appearing in 
each. We replace the t integration with an integration over 4 variables i; (i = 1, ... ,4), replacing t with ti in each 
term in the denominator, and adding in delta function constraints to make all the ti equal. This is just a rewriting of 
the original integral. We then use the integral representation of the delta functions to obtain: 



P(u,x) 



2tt 



du>idu)2(lu)3 exp 



r(2/i5) 2 r(i/i5)V 17 /i5 Uc i6/i5 

"13/15, ,"13/15, I-W15,, ,_|-14/15 Q(w + _ ^(wa _ ^2)6(^3) (19) 



\<JJ + OJ\\ 



\LU 2 - LUl\ 



W 3 - LU 2 \ 



N3 



Because the integrals over ojj are all convergent in the above expression, we see that P(to, x) is continuous as a function 
of x. We write 



With this notation we obtain 

ne 



U(t)) = 



5r(2/5) 



P(lO,x) = P(w,0)H(LJx/Vn,V n /Vc) 



3/5 sgn(Luj)v n - 4 / 15 v c - 2 / 15 C jk T j T k *H{u J {x j - x k )/v n ,v c /v n ) 



(20) 
(21) 



IV. FINITE TEMPERATURE 



The analysis of finite temperature effects again proceeds much in the same way as in [5j . All the nontrivial topological 
information is still encoded in the Cjk- The only difference between the zero temperature and finite temperature case 
arises from the fact that the 2-point functions we now wish to use are analytic continuations of 2-point functions on 
a cylinder. The replacement that needs to be made is^: 



[5 + i{t±x/v n )] 



2/15 



ttT 



sin[7rT((5 + i{t±x/v n ))] 



2/15 



(22) 



for the <7i correlator, and similarly for the correlator of the exponentials of the charge bosons, except with v n replaced 
by v c . We obtain the following modified P: 



P(w,x,T)= / (dt) 



sin[7rT(<5 + i(t + x/v n ))] 



ttT 

sin[7rT(5 + i(t + x/v c ))] 



-2/15 



sin[7rT(i5 + i(t — x/v n ))] 



ttT 



-2/15 



ttT 



-1/15 r • 



sm[TrT(S + i(t-x/v c ))] 



-1 -1/15 



ttT 



(23) 



As opposed to the case of [5]. this integral is not exactly doable for arbitrary values of v n ,v c , the neutral and charge 
mode velocities. However, we can extract the behavior of the current at large values of the temperature. To do this, 
simply note that the integrand above is basically a product of inverse powers of sinh, which at T large compared to 
hui^ sc and hoj" sc become decaying exponentials, suppressed everywhere except at t — ±x/v n ,±x/v c . Noting that the 
integral is everywhere convergent, we see that most of its contributions come from those 4 points. Plugging in, we see 
that we have an exponential suppression of 



exp 



-2t:T\xi — X2I 



2/15 1/15 



(24) 



of the interference term at high temperatures T 3> Tiuj^ sc1 Tilo™ bc . Comparing this to the result of Q who obtain 



exp [—2irgT\xi — x 2 \ 
we see that we have an "effective" g of 2/(15v n ) + l/(15i; c ). 



(25) 
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V. DETERMINING THE C jk 

We now come to the crux of the problem, which is to find the constants Cjk- This basically comes down to 
determining the conformal blocks for the Oj correlators, and here the topology comes into play. Working in the 
context of the v — 5/2 state and a single point contact, Q argued that for a given quasihole tunneling process, 
the two cr fields that create and destroy the quasihole should fuse to the identity. Arguing along these lines, one 
can uniquely determine the desired conformal block. It is here that one uses the connection between the bulk and 
the boundary: the tunneling topology in the bulk determines the fusion channels of the tunneling operators in the 
boundary field theory. To give an actual argument for this connection would require a detailed examination of bulk 
tunneling, which has not been done. We shall, as in Q, adopt it as part of the definition of the tunneling operator, 
and adapt it to our case of a double point contact. 
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FIG. 1: Tunneling geometry with a test quasihole between the two point contacts. The test quasihole is created and annihilated 
on the same side of the hall bar 

We wish to probe the effect of test quasiholes, residing on anti-dots in the region between the two point contacts, on 
the tunneling current. The main conceptual point we wish to stress is that to do the computation correctly, one must 
incorporate the creation and annihilation of these test quasiholes in the far past and future. That is, we must include 
two more aj operators - we will take them to be on the same side of the hall bar, creating the quasihole in the far past 
and annihilating it in the far future (see fig l[T))). For simplicity, we will analyze just two situations: one without any 
test quasiholes, and one with a single test quasihole, and for easy comparison we will actually create a test quasihole 
in both situations, but keep it far away from the region of interest in one case (so that it has trivial braiding). For 
definiteness, we will for now analyze the following correlator, one of two that contribute to the interference term: 

(a 2 (xo 1 t )a 2 (xi 1 t)ai{D - xi,t)a 2 (D - x 2 ,t')a 1 {x2,t')a 1 {x , -t )) (26) 

Here the test quasihole is created at (xq, —to), and annihilated at (xo,to). to and D are to be taken to infinity. We 
find it convenient to send D to infinity first, so the thing splits into a 4 point and a 2 point function by cluster 
decomposition. Let's see what happens when to goes to infinity. Then <J2{xi,t) and a\{x2 1 t ! ) can be considered to 
be very close (in comparison with (72(2:0, to) and <Ji(xq, —to)), so we can use the operator product expansion on them. 
Clearly the leading contribution comes from the identity, so that's the only one we need to retain in the limit. The 
4-point function is then just the 2-point function of 02(2:1, t) and 01(2:2, t') times an overall constant having to do 
with the preparation and annihilation of the test quasihole. 

From the above paragraph, we see that it is very convenient to work in the basis where the fusion channels of the 
pairs (02 (a;o> to)) 0i (#0 j — to))i (02(^11 t), o\ (#2, t')), and (01 (D — x±, t), 02 (-D — x — 2, f')) are known, because only in 
the case when all those fusion channels are the identity do we get a nonzero contribution. But the conformal block 
we're interested in is specified in a different basis. To express it in the new basis, we proceed as in Q and project the 
worldlines of the quasiholes (including the test quasihole) to two dimensions; the linear transformation necessary to 
enact the basis change is the matrix representing the resulting braid. 

The conformal block we start with is: 



FIG. 2: Projection of the quasihole worldlines in the case of a test quasihole braiding nontrivially with the others 




FIG. 3: Projection of the quasihole worldlines, together with a test quasihole braiding trivially 



(X2,f) (D-X2,f) (Xl,t) (X2,f) (xo,—to) (X0,t()) 
I | | I | | I | | J 

We want to expand it in the new basis above, and in particular get the coefficient of 



(27) 



(X2,f) (Xi,t) 



(D-x 2 ,f) (x 2 ,t>) 
I I I 



(x ,-t ) (x 0) t ) 

I I I / 



(28) 

Depending on whether or not the test quasihole resides in the region between the two point contacts, we have 
two possible braiding histories, shown in figures @ and ([3]). Because, due to the considerations outlined below, we 
only care about the coefficients up to phase, we only need to determine the matrix for the linear transformation 
between the two bases up to phase. Thus, even though a priori we would have to deal with the Moore-Seiberg braid 
matrix for the parafermions (and the associated U(l) theory), we really only need the braid matrix in the associated 
quantum group picture. Here, both o\ and a 2 correspond to the 2-dimensional spin 1/2 representation of U q (sh), 
q — exp(27ri/5). Labeling conformal blocks using the Bratteli diagram notation, we see that we start off with the 
block represented by the path [(0, 0), (1, 1), (2, 0), (3, 1), (4, 0), (5, 1), (6, 0)], which for convenience we will designate via 
shortened notation by [0101010]. We see that when the braids in ([2]) and |3]) act, everything happens in the subspace 
spanned by [OlxlylO], x,y = 0,2, and each transposition acts either only on x or only on y (they commute). The 
relevant unitary matrix enacting the transposition is [||: 



Ui. 



— 7H/10 / -27ri/5 



i+1 



(29) 



where i = 2,4. The 2-element basis in which we write the matrix t/2,3 is x = 2,0, and the one corresponding to U4.5 
is y = 2,0. Note also that, using some identities involving the golden mean, we have: 



We have 



W5 ^ e 4W5 i^ T (j + 2) 

i^T(j + 2) -r 4 "/ 5 



(30) 



([0101010]|[/ 2 , 3 |[0101010]) 
([0101010]|[/ 2: 3^ 2 4, 5 | [0101010]) = 



(31) 
(32) 
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Let us now use this to obtain the Cjk- First, let's see what we can surmise on general grounds. By tuning T2 to 
and using the fact that (j(t)) is real, we find that Cn is real. Also, by symmetry, Cn = C22 = C. Again because the 
current is real we see that C21 — C12* , Cyi = \C12\e . The phase 9 can in principle be determined, but we ignored it 
in the calculations above. This is because the effect of changing 9 is the same as that of changing the phases of Fi 
and T2, which are hard to determine anyway and are sensitive to the Aharonov-Bohm effect. 

The results in (f3"Tj) and (|32[) show that when a quasihole is present, | C12 1 is diminished by a factor of 1/r 2 . Similar 
topological calculations (with only the spatial coordinates of the points involved changing) show that the other 
constants, C\\ and C12, are always equal, and equal to the | C 12 1 in the case when there is no test quasihole in the 
region between the point contacts. Thus we have 

(i(*))noqh. = A\LOj\- 3/5 S gn(uj)v n - 4 / 15 v c - 2 / 15 (Vxl 2 + |r 2 | 2 + {T^* + T 1 *T 2 )H{uj\x l - x 2 \/v n ,v c /v n j) (33) 



(?'(*)> q h. = A\uj\- 3/5 sgn(u;j)v n -V 15 v c - 2 / 15 (jri| 2 + |r 2 | 2 + ^(T,T 2 * + r^r^M^ - x 2 \/v n ,v c /v n )J (34) 

where A is an unimportant constant. 

We thus observe that even without looking at phases, nonabelian statistics have a robust effect of changing the 
magnitude of the coefficient of the interference term in the tunneling current. This is in line with earlier arguments 
of [Tj for the differential conductivity. 
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